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Abstract 

O 

Motivated by a conjecture put forward by Abramowicz and Bajtlik we recon- 
sider the twin paradox in static spacetimes. According to a well known theorem 
in Lorentzian geometry the longest timelike worldline between two given points 
is the unique geodesic line without points conjugate to the initial point on the 
£>. | segment joining the two points. We calculate the proper times for static twins, for 

twins moving on a circular orbit (if it is a geodesic) around a centre of symmetry 
and for twins travelling on outgoing and ingoing radial timelike geodesies. We 






show that the twins on the radial geodesic worldlines are always the oldest ones 
^ ' and we explicitly find the the conjugate points (if they exist) outside the relevant 

segments. As it is of its own mathematical interest, we find general Jacobi vector 
fields on the geodesic lines under consideration. In the first part of the work we 
investigate Schwarzschild geometry. PACS number: 04.50.Kd 
Keywords: twin paradox, static spacetimes, Jacobi fields, conjugate points. 
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1 Introduction 

The twin paradox is well understood in Minkowski spacetime whereas in curved space- 
times less is known about generic features of the twins following various worldlines and 
recent works studying various special cases are quite numerous (see e. g. [Tj, [2J, [3J, [I], 
[5J, [B] and references therein). Abramowicz, Bajtlik and Kluzniak investigated a couple 
of static spacetimes where the absolute rest is denned by the unique timelike Killing 
vector field and a stationary axially symmetric spacetime where the standard of rest is 



determined by the rest frame of the zero angular momentum observers (ZAMO) and 
calculated the proper times for the twins orbiting on the same circular orbit with dif- 
ferent constant linear velocities. From the case of these circular and stationary motions 
they have drawn a conclusion which being generalized to other worldlines is actually 
a conjecture: 'in all situations in which the absolute motion may be defined in terms 
of some invariant global properties of the spacetime, the twin who moves faster with 
respect to the global standard of rest is younger at the reunion, irrespectively to twins' 
accelerations' [7] , [8] . It is the purpose of the present work to show that in general the 
conjecture is false. 

We first make a comment on the concept of rest. Geometrically an absolute standard 
of rest is provided by the timelike Killing vector field (if it exists and is unique), one may, 
however, consider a universal standard of rest based on the synchronous (i.e. comoving, 
i.e. normal Gauss) frame of reference which is freely falling in the gravitational field. 
The comoving frame may cover only a part of the spacetime (as it develops, besides some 
special cases, coordinate singularities) and the distances between freely falling bodies 
change in time, nevertheless to some extent the bodies may be treated as being at rest. 
Accordingly, if radially freely falling bodies in Schwarzschild spacetime are supposed to 
be at rest, then a static Killing observer and an observer on a circular geodesic (i.e. an 
observer orbiting around the central body) are moving with respect to them and, as is 
shown below, their proper times are shorter than those at rest. This case supports the 
Abramowicz's, Bajtlik and Kluzniak (ABK) conjecture. However, this notion of rest is 
not unambiguous. In fact, the observers on circular orbits are freely falling too and if 
they are regarded as staying at rest then the radially falling observers and the static 
Killing ones are in respective motion and their proper times are larger than those of the 
motionless (i.e. orbiting on circles) observers. We therefore assume that the state of rest 
is determined by the unique Killing field. 

The problem of which twin travelling between two given spacetime points is older is 
of purely geometrical nature and the answer to the question of which worldline makes a 
twin the oldest one (assuming a large number of twins following various worldlines with 
the same endpoints) is implicitly well known as it may be found in advanced mono- 
graphs. Imprecise answer is even contained in an intermediate level textbook [9]: 'for 
a really reliable answer one has to know how to deal with accelerated systems; here 
general relativity is to be asked, and answer is: yes, travelling (deviating from geodesic 
motion) keeps you younger'. 

A rigorous formulation and solution of the problem is following. Let U be a convex 
normal neighbourhood of a point p in any Lorentzian spacetime. Then the following 
proposition holds: 
Proposition 1 (proposition 4.5.3 in [TO] ) 



In any convex normal neighbourhood U, if p and q can be joined by a timelike curve 
then the unique timelike geodesic connecting them has length strictly greater than that 
of any other piecewise smooth timelike curve between the points. 

Schwarzschild spacetime is misleading in this context. Consider two twins: A is at 
rest at a point r = r$ > 3M (in the standard Schwarzschild coordinates) and B moves 
on a timelike geodesic line with a circular orbit in the space at r = r . Initially the 
twins are at a spacetime point Pq (t = t , r = r ), then B flies away, makes a full circle 
around the centre and meets A again at Pi (t = ti, r = ro). A simple calculation shows 
that the proper time sa (i.e. the worldline length) of A between P$ and Pi is longer 
than the geodesic proper time sb of B, 
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where v is B's velocity with respect to the static A [7], [8]. This case gave rise to the 
conclusion that the twin who moves faster (and without acceleration) is younger than 
the static twin (subject to acceleration). However this simple example merely shows 
that Pi does not belong to a convex normal neighbourhood of Pq. In fact, there is at 
least one timelike geodesic line joining Po and Pi whose length is larger than that of B. 

Here the key notion is that of conjugate points and for the reader's convenience we 
briefly recall the relevant facts. Let 7 be a timelike geodesic with tangent unit u a . Any 
vector field Z^(s) which is a solution to the geodesic deviation equation on 7 is called a 
Jacobi field on 7. A pair of points p and q on 7 are said to be conjugate if there exists a 
Jacobi field Z^ on 7 which is not identically zero and Z^(p) = Z^(q) = ([10J chapter 4). 
Points p and q are conjugate if infinitesimally nearby geodesic intersects 7 at both p and 
q. Existence of conjugate points on a geodesic shows that it is neither the unique geodesic 
nor the longest curve connecting its endpoints. More precisely, if a timelike geodesic 7 
joining points p\ and P2 has a point q conjugate to p\ belonging to the segment P1P2, then 
there exists a nearby timelike curve r (not necessarily a geodesic one) with endpoints p\ 
and pi, which is longer than 7, s{j) > s(j). Conversely, if there is no conjugate points 
on a timelike geodesic, then it is the longest curve between its endpoints, independently 
of whether the endpoints lie in some convex normal neighbourhood, or not. It is stated 
in 

Proposition 2 (proposition 4.5.8 in [TO] ) 

A timelike geodesic 7 has the maximal length from p\ to P2 if and only if there is no 
point conjugate to p\ on the segment p\p<z- In some spacetimes, e.g. anti-de Sitter space, 
two given points may be connected by various timelike curves, but none of these attains 
maximal length. The longest curve always exists in globally hyperbolic spacetimes. 
Proposition 3 (proposition 9.4.4 in [TO] ) 
Let (M, g a p) be a globally hyperbolic spacetime. Let p\ and pi be connected by timelike 



curves. Then there exists a timelike geodesic 7 from pi to p<i having the maximal length. 

The existence of conjugate points is determined by 

Proposition 4 (proposition 4.4.2 in [TO] ) 

If R a p u a u 13 > on a timelike geodesic 7 and if the tidal force R^ au p u a u 13 7^ at some 

point po on 7, there will be conjugate points p and q on 7, providing that the geodesic 

can be extended sufficiently far. 

Returning to the problem of the rate of ageing of the twins one sees that neither the 
velocity (in stationary spacetimes) nor absence of acceleration (on geodesic worldlines) 
alone is sufficient to establish which twin will be older at the reunion. In a curved 
spacetime there is a great variety of motions and for a given pair of worldlines there is 
no generic criterion (if none of them is a geodesic free of conjugate points) stating which 
of these is longer; one must explicitly compute their lengths. Therefore besides studying 
case by case one can only determine the longest worldlines joining given two points in 
some spacetimes. This is the heart of the present work. 

For a given pair of spacetime points one applies physical arguments to find out a 
timelike geodesic connecting them and seeks for conjugate points on it. In physically 
interesting spacetimes the strong energy condition holds implying R a @ u a u 13 > and in 
most cases the tidal forces do not vanish, hence Proposition 4 indicates that the geodesic 
under consideration contains somewhere conjugate points. It is then crucial to establish 
whether the conjugate points belong to the relevant segment of the geodesic. 
In this work we explicitly analytically compute location of conjugate points on physically 
distinguished timelike geodesic lines in a number of static spacetimes and for comparison 
in Robert son- Walker spatially flat universe. To establish that the given geodesic is not 
the maximal one it is sufficient to find out a particular Jacobi vector field on it which 
reveals conjugate points within the relevant segment. A more general problem is finding 
out a generic Jacobi field on the given geodesic curve and there is an opinion (private 
communication) of a few experts in the field that this work is worth doing it as it is 
of mathematical interest. We get general solutions of the geodesic deviation equation 
for worldlines followed by the twins considered in the problem of their ageing rates. In 
the first part of the work we discuss the most convenient form of the geodesic deviation 
equation for the Jacobi field and then deal with Schwarzschild spacetime. 

2 Timelike geodesies in static spacetimes and Jacobi 
fields 

We consider static spacetimes since there the orbits of the timelike Killing vector field 
are orthogonal to hypersurfaces of simultaneity and the notion of observers at rest is 
unambiguous. The case of standard of rest with respect to the ZAMO frames shows that 



stationary non-static space-times require a separate treatment [5] . In a static spacetime 
there always exists a coordinate chart in which the metric takes the form 

ds 2 = g o{x k ) dt 2 + gij{x k ) dx i dx j (2) 

and the timelike Killing field is K a = <J«, K a K a > C0, a, /3 = 0, 1, 2, 3 and z, j, fc = 1, 2, 3. 
An observer is at rest if his worldline coincides with one of the orbits of the Killing field, 

-| j<y 

then his velocity is u a = g 00 5q. From the geodesic equation (D/ds)u a = 0, where 
D /ds is the absolute derivative along a curve with respect to its arc length, one imme- 
diately sees that a static observer worldline is a geodesic if and only if g 00 — 1- This is 
a very special case and none of the physically interesting static spacetimes satisfies this 
condition. We shall always denote by A the nongeodesic observer at rest. 

We assume that the spatial coordinates (x l ) are so chosen that (d/dx^goo > 0. The 
coordinate x l varies in the range — oo<a<x 1 <6<+oo and x l tending to b corre- 
sponds to spatial infinity (if it exists). Due to these properties we dub x l into 'radial 
coordinate'. We then seek for two kinds of worldlines as possible geodesies. Firstly, 
'circular orbits' with x l = const ^ which are spatially closed and secondly 'radial 
worldlines' with x 2 and x 3 constant and x 1 varying in an subinterval ao < x 1 < bo of 
(a, b). In the spacetimes studied in this work the radial curves are timelike geodesies. 
In the physical interpretation we shall assume that the twin B follows a circular world- 
line and the twin C moves on a radial geodesic curvcl The twin A stays at x 1 = a and 
in a coordinate time interval At his proper time is sa = y/goo~At. The twin C's radial 
worldline emanates outwards from x 1 = ao, reaches a maximal height x 1 = bo, turns 
inwards and comes back to x l = a to meet again A and possibly B. The corresponding 
proper time Sc of the twin C may be greater than sa] this is possible although the 
expression g^ dx l dx^ < for dx l ^ if goo grows sufficiently quickly in the interval 
ao < x 1 < bo. 

On a given radial or circular geodesic 7 the conjugate points can be determined 
directly from the definition: by finding a generic solution Z^(s) to the geodesic deviation 
equation 

ds 

where u a = dx a /ds, u a u a = 1, is the tangent vector field to the geodesic. The Jacobi 
field is orthogonal to the velocity vector, Z a u a = 0. Due to the presence of the absolute 
derivatives, equations (3) are very complicated ODEs of second order and to simplify 
them one applies the following approach ([10J chapter 4 and [11] chapter 9). On 7 one 



2 Z» = R\ M u a uPZ\ (3) 



^or obvious reasons we choose the metric signature (H ). 

2 We use here the traditional term 'twins' instead of the correct 'triplets'. 



introduces an orthonormal triad of spacelike vector fields eg(s), a = 1,2,3, which are 
orthogonal to m m , 

<%e bli = -6 ab , e> M = 0. (4) 

Z^ are components of the Jacobi field in the coordinate basis, Z = Z^d^. Since Z a u a = 
0, the field may be decomposed in the triad basis, Z = ^ a Z a e a and when the triad 
fields are kept fixed the functions Z a (s) are three scalar fields on 7. Equations (3) get 
simplified when eg are parallelly transported along 7, (D/ds)e% = 0. The components 
Z^ are expressed in terms of the scalars as 

^ = X>aC (5) 

o=l 

Inserting (5) into (3) and multiplying the result by eg one arrives at the following linear 
equations for the three scalars Z a , 

^Z a = -eg R mM u a U? J2 e< Z b . (6) 

aS 6=1 

At an initial point po on 7 one gives initial values eg(po) satisfying (4), then the triad 
fields are uniquely determined by the parallel propagation. One seeks for a general 
solution to (6); for each a it depends on two integration constants and one of them is 
a multiplicative factor. In the search for conjugate points one then imposes the initial 
condition Z a (p ) = and gets the desired Jacobi field. If one of the scalars, say Z\, 
vanishes at a point q on 7, Zi(q) = 0, one can identically set Z 2 = Z 3 = and the 
Jacobi field Z fJ, (s) = Zi(s) e^ (s) is zero at q; the points p and q are conjugate on 7. 

In the next section we calculate and compare the lengths of a circular and a radial 
geodesic worldlines, determine Jacobi fields on them and their conjugate points in the 
case of the static spherically symmetric vacuum spacetime. 

3 Schwarzschild spacetime 

In this spacetime (in standard coordinates t, r, 6, <ft) geodesic circular orbits exist only 
for r = r > 3M, where M is the energy of the spacetimqj. All the motions take place 
in the two-surface 9 = ir/2. 

The twin A stays at rest at r = vq, 9 = tt/2 and = <f>Q. The twin B moves on a 
geodesic circular orbit r = r$ with a constant linear velocity vb with respect to A. All 
the parameters of the B's geodesic are determined by r (and M) and it is parameterized 



3 Usually we put c = G = 1; from time to time we shall explicitly display c for clarity. 



by its length, 

t - to =(-l^f s , ,_, .IpiLV /2 s (7) 

\r - 6MJ r Q \r - 3M J 

and (vb/c) 2 = M/(r — 2M); t is the moment when A, B and C are together at the 
same place. The length of the geodesic B after making one full circle is determined from 

(7) as — 0o = 27r and is 

/r -3M\ 1/2 
s B = 2nr -5-—) (8) 



M 
and the corresponding coordinate time interval is 

A( = 2 "'»S) I/2 - (9) 

By assumption at the point P^it = t ,r = r , = O ) the worldlines of the three twins 
intersect. The twin C moves on a radial geodesic emanating from Pq with an initial 
radial velocity (dr/ds)(t ) = u > directed outwards at 6 = tt/2 and = 0o- At £ = tu 
he reaches a maximum height r = tm where his radial velocity decreases to zero and he 
falls down along the ingoing radial geodesic until at t — t\ he returns to r = r . Since 
the spacetime is static, the coordinate and proper time of the flight upwards are equal 
to these of falling down, i.e. t\ — tu = tu — to- 

By assumption the three twins meet again at P\{t — ti,r — r , = 0o), i.e. C falls back 
to the initial place when B makes a full circle, 

h-t = 2{t M - t ) = At. (10) 

The proper time sa of the static twin from P to P± is 

/ 2M\ 1/2 /r n \ 1/2 

* = (i--) *=^GiH (11) 

and according to (1) is larger than sb [I], 13, what shows that the circular geodesic B 
contains conjugate points between P and P\. 

One then conjectures that the longest curve connecting Pq and Pi is the radial geodesic 
C. To confirm the conjecture one computes the C's length and then shows that conjugate 
points, if exist on C, are located beyond the segment PqP\. As is well known it is 
convenient to parameterize a radial geodesic in this spacetime by an angular parameter 

V, 

r{ri) = r M cos 2 77 = -r M (cos2r/ + 1), (12) 

where —a/2 < rj < tt/2 and tm is the largest radial distance from the centre. By 
definition the initial point is r(— a/2) = r , then 

2 r 
cos a = 1 and < a < it. (13) 



The return point is r(a/2) = r , the event horizon r = 2M is at i] = rjn where cos 2 t]h = 
2M/tm and r(ir/2) = 0. On the radial geodesic all geometrical quantities are determined 
by tq and the initial velocity u, or equivalently by the integral of energy. The twin C 
with mass m and 4-momentum p a has conserved energy E determined by the timelike 
Killing field, which is normalized to unity at the spatial infinity, K a = Sq, hence 
E/c = K a p a = p = (1 - 2M/(c 2 r))mcdx°/ds > 0. 
We introduce dimensionless constant energy per unit rest energy, k = E/(mc 2 ), then 

where x° = t. The geodesic equation for the radial coordinate is replaced by the universal 
integral of motion u a u a = 1 which in the present case reads 

,o 2M 
= k 2 + 1. (15) 




At the starting point one has 



« 2 = fc 2 + — -1 (16) 



and tm is determined from k 2 — 1 + — = 0, hence 

rM = ~\ Z2 and < /c < 1. (17) 

Since stable circular orbits require Tq > 3M and tm > r o one gets a lower limit k 2 > 1/3. 
The tangent vector u a has, from (14) and (15), the nonvanishing components 

dt k cos 2 r) dr 



— = -Vl- k 2 tanr^; (18) 

ds cos z 7] — 1 + k A ds 

the minus sign ensures that for —a/2 < rj < the velocity is directed outwards and for 
i] > it is directed inwards. To obtain relations t = t(rj) and s = s(i]), one replaces (18) 
by derivatives with respect to rj. The derivative ds/dr] is derived from the formula for 
ds 2 expressed as a function of rj; one gets 

= ^cos 4 r? (19) 

and employing that ds/dr] is always positive one finds 

/„3 \V2 r 




s(rj) 



1 M 

2M 



1 . . . o; 

-(sm2?7 H-sinaJ +77 + — 



(20) 



This immediately yields the length of C from Pq to Pi 

/ j& \i/2 r 



sc 



a 



• M 

2M 



2r r M (2r 

a / 1 + arccos 

Tm V r o \r M 



(21) 



As concerns t{rj) one has dt/drj = (dt/ds)(ds/drj) and substituting from (18) and (19) 
one gets 



t(v) ~ to 



r 



1/2 



M 



1 M 

2M 



2M 

1/2 



(cos2r7 + l) 2 
cos2t? + 2A; 2 - 1 



c/r; 



a 1 
[(3 - 2A; 2 )(r7 + - ) + -(sin 2?/ + sin a)] + 



+ 2M In 



6 + tan?/ 
6 — tan r\ 



2Mln 



6 — tana/2' 
b + tana/2 



(22) 



wherqfl 6 = fe(l — fc 2 ) _1//2 . Notice that at the event horizon sin tjh = +k, then t&nrjH = b 
and t — y +oo as it should be. Accordingly, the coordinate time interval of the flight 
back and forth is 



1/2 



t(a/2) - t 



' M 

2M 



[(3 - 2A; 2 )a + sin a] + AM In 



b + 



' Tm 

- r 



1)1/2' 



'TM. _ 1)1/2 



(23) 



Here 



V - I r JL - 1 



1 



1 



2M 



>0, 



r / 1 — A; 2 V r 

hence the denominator of the fraction under the logarithm is always positive. The in- 
terval (23) must be equal to At given in (9) and this gives a transcendental algebraic 
equation for the energy k or equivalently the initial radial velocity u in terms of the pa- 
rameter ro- Introducing for convenience a parameter R = ro/(2M) > 3/2 and replacing 
the unknown fc(ro) by x = tq/tm = R (1 — k 2 ) = 1 — u 2 R one finally gets an equation 
for x(r ) (or x(R)), 



2V2nR 3/2 



x 



1/2 /Rx 3/2 

/? / V x ) 



2x 
~R 



a 



2V. 



x — x A 



+ 



+ 2 In 



y/R — X + a/1 — X^ 



Vr 



X 



VT 



x . 



(24) 



where, according to (13), cos a = 2x — \. It may be shown that for each value of R > 3/2 
there is only one solution of the equation in the interval < x < 1. Inserting this value 



4 Clearly formulae (20) and (22) have already been known in various versions, we rederive them here 
for the reader's convenience. 



Table 1: The ratio sq/sa for 5 values of x(R). 



R 


X 


sc/sa 


1,5 


0,545246 


1,19878 


1,6 


0,534476 


1,17554 


2 


0,508559 


1,12061 


4 


0,472995 


1,04849 


100 


0,44932 


1,00162 



of x(R) into (21) one gets the ratio of the twins' proper times, 

S C 1 (, l Y 1/2 -3/2 [ I 2 , 1 

— = = 1 x ' vi — x z H — arcci 

sa v^vr V i?/ I 2 



Vx — x 2 H — arccos(2a; — 1) 



(25) 



Equation (24) was numerically solved for five values of R consecutively growing by 
0,1 from 1,5 to 2 and for all integers from R = 2 to 100. The value of x(R) is a 
monotonic slowly diminishing function of R (we did not attempt to fit it by an analytic 
expression) and decreases from 0,545246 for R = 1, 5 to 0,44932 for R = 100. It is hard 
to analytically establish from eq. (25) which proper time is larger and we compute their 
numerical values. In table 1 we give the ratio sc/sa for some chosen values of x(R). As 
it should be expected, the farther the point r is from the event horizon the smaller the 
ratio Sc/sa is and always the radial geodesic C is longer than the non-geodesic curve 
A. That s c > sa was previously found in a special case in [T2"] . 

Both the geodesic curves B and C satisfy the conditions of Proposition 4, hence 
there do exist conjugate points and to identify them we now find Jacobi fields on these 
worldlines. 



3.1 The timelike radial geodesic C 

An alternative expression to formula (18) for the tangent vector to C in terms of the 

coordinate r is 

kr (2M -\ 1/2 



u 



r-2M' 



-l + k x 



.0,0 



(26) 



where e = +1 on the outwards directed segment of C and e = — 1 on the inwards 
directed piece. The parallelly propagated and orthogonal to C orthonormal spacelike 
triad is chosen as 



( i 



er 



2M , 9 \ 1/2 
1 + k 2 



r -2M \ r 
_(1_£;Y /2 sinr ? cosr ? 



COS 2 T] — 1 + k 2 
10 



,fc,0,0 
,fc,0,0 



(27) 



c 2 



0,0,-, 

r 



e 3 



0,0,0, 



(28) 



with e as in (26). The geodesic deviation equation for the scalars Z a reduces in the 
present case to 

*- 2M„ (2g) 

(30) 

(31) 

Simplicity of these equations is deceptive since their left hand sides contain derivatives 
with respect to the arc length instead of r or 77. One replaces these second derivatives 
by derivatives with respect to 77, 



d^^ 1 ' 


r 3 ~ 1 ' 


d 2 
d? Z2 = 


M z 


d 2 


JL Zz . 



dZa 

ds 



P 



dZa 

cos 2i] + 1 drj 



cPZ a 

ds 2 



jr 



(cos2?? + 1) 



d 2 Z a 2sin2r] dZ a 
dr) 2 cos 2i] + 1 drj 



where p = {r\j/2M) x l 2 . Then the Jacobi equations read 

d 2 Z 1 2sin2r] dZ ± 8 



+ 



(32) 
(33) 

(34) 

(35) 
drf ' cos 2r/ + 1 dr/ ' cos 277 + 1 

and for Z3 one has an equation identical to that for Z2. To solve eq. (34) one introduces 
a change of the independent variable, x = cos 2r) (not to be confused with the parameter 
r o/ r M introduced earlier), then u — (x + l)/2 and substitutes Zi(u) = 4u 2 G(u). Then 
(34) reduces to a special case of Gauss hypergeometric equation which is solved in 
elementary functions. After some manipulations one gets a generic solution for Zi, 



drf cos 2?7 + 1 dr] cos 2^+1 



d 2 Z, 



+ 



2 sin 2r] dZ-i 



+ 



£i = 0, 







Z\ (77) = C\ tan 77 + Ci (3?7 tan 77 n tan 77 — cos 2 77 + 3) . 



(36) 



In the case of eq. (35) one performs the same operations as previously and finally 
substitutes G(u) = - Y(u), this yields an equation 



, d 2 Y ( 3 



dY _ 

du 



(37) 



11 



which is also solved in elementary functions. The general solution is 

Z 2 (r]) = C sin 2t? + C"(cos 2r) + 1) 



(38) 



and an analogous formula holds for Z3. A general Jacobi field on a radial geodesic was 
earlier found by Fuchs [13] . He used a different basis triad whose explicit form in terms 
of r or 1] was not given, hence it is not easy to compare his solutions to ours. 

We now seek for conjugate points on C. The special Jacobi scalars which vanish at 
Po(t] = —a/2) are then 



Zi(v) = c x 



3 4 cos a — cos a + 5 \ 2 

-a -\ h 37? tan 77 — cos 77 + 3 

,2 2 sin a j 



z 2 {v) = c 2 



sin 277 + 



ana 



(cos2t7 + r 



(39) 



(40) 



cos a + 1 

and Z 3 is given by the same formula as that for Z 2 with C 2 replaced by C3; C\, C 2 and 
C3 are arbitrary constants. Functions Z 2 and Z3 vanish at 77 = n/2, i.e. at r = 0. In 
this sense the curvature singularity is a point conjugate on the radial geodesic to any 
point P . For —a/2 < r\ < ix j2 the equation Z 2 {r\) = is reduced to 



tan r] 



sin a 



cos a + 1 



?"0 



and its solution belongs to the interval —tt/2 < rj < where tangens is monotonic, hence 
the only solution is 77 = —a/2. In sum, the Jacobi field Z^ = Z 2 e 2 + ^363 is different 
from zero on the radial geodesic segment from any P to the event horizon. 
As concerns Z\{r\) one easily shows that the ratio Z\/C\ > 2 in the interval < 77 < n/2 
and possible zeros may be only on the outgoing segment of the geodesic, —a/2 < 77 < 0. 
One views Z\ as a function of two variables, Zi(j],a) with a G (0,7r). For any fixed 
value of 77 one has a G (—27/, it) and in this interval Zi(rj, a) grows since 



dZ x 
da 



cot 2 a/2 cos 2 a/2 tan(— 77) > 0. 



From Zi(—a/2,a) = one infers that Zi/C\ is positive for —a/2 < rj < 0, thence the 
function has no zeros for —a/2 < 77 < %/2. 

This completes the proof that a timelike radial geodesic consisting of an outgoing 
segment and an ingoing one contains no conjugate points and is the longest curve joining 
points Pq and Pi outside the event horizon. Clearly the same statement concerning 
non-existence of conjugate points applies to any timelike radial geodesic which is purely 
outgoing or ingoing. 
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3.2 The timelike circular geodesic B 



It is rather surprising that solving the equations for Jacobi fields is in this case much 
harder than for a radial geodesic. One immediately gets the vector u a tangent to the 
geodesic B from its parametric form (7). A triad of spacelike orthonormal vector vector 
fields on B satisfying (4) and parallelly transported along it may be chosen as 



e i 



Mrn 



1/2 



(r -2M)(r -3M) 
-1 /r -2M\ 1/2 



sin qs, 



r \r — 3M 



sings 



p' 1 

1 2 



1 
r 



0,0,-, 



r 



2M\ 1/2 



r 



cos qs, 0, 



(41) 



(42) 



p' 1 

1 3 



Mr n 



(r -2M)(r -3M) 
r -2M\ 1/2 



1/2 /r -2M\ 1/2 

cos qs, I I sings,0, 



r \r — 3M 
where q 2 = M/r^. Employing 



cosgs 



u 



(r - 3M)" 1 / 2 (rj^ + ^5? 



and denoting 



/3 



r -2M 



r -3M 
one arrives at the following form of eqs. (6) for the scalars Z a (s), 

d 2 
-— ■Z 1 = q 2 [(j3 cos 2 qs — 1) Z 1 + (3 Z% sings cos qs], 



ds 1 



M 



r 2 (r - 3M) 



ds 2 



Zz = q 2 \ft Z\ sin qs cos qs + ((3 sin 2 qs — 1) Z3] . 



(43) 



(44) 



(45) 
(46) 
(47) 



It will be shown that Jacobi fields spanned on the basis vectors e^ and e^ do not 
give rise to a point conjugate to Pq on the segment PqP\ and in this sense they are 
irrelevant to the twin paradox. We therefore postpone solving equations (45) and (47) 
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to the Appendix and in the current subsection we discuss consequences of the equation 
for Z 2 - 

The general solution for eq. (46) reads 

Z 2 = C sinp s + C" cospos (48) 

where 

Pl = rg(r - 3M) (49) 

and C and C" are integration constants. The special Jacobi field proportional to e 2 and 
vanishing at Pq{s = 0) is simply Z^ = C 5 2 sinpos. The field generates infinite number 
of points conjugate to Pq which are located at 

- = "' (t)"^ < 50 > 

n = 1,2, . . .. The nearest to Pq conjugate point Q is at sq = sb/2, that is half way 
between Pq and the endpoint P±. Notice that Pi is also conjugate to both Pq and to Q. 

One conjectures that analogous features occur in Reissner-Nordstrom spacetime, 
this will be shown in a forthcoming paper. 

4 Summary 

In this work we give explicit forms of generic Jacobi fields and conjugate points gener- 
ated by these fields on physically distinguished timelike geodesic curves in Schwarzschild 
spacetime. These are circular orbits and radial curves going back and forth. Applying 
the theorem that the longest timelike curve joining two given points is a geodesic free 
of points conjugate to its endpoints and lying between the ends, we show that in this 
spacetime the longest worldline is any radial geodesic. This outcome may be physically 
interpreted in terms of the famous twin paradox in the spacetime. Considering a number 
of twins following different worldlines with common endpoints one finds that contrary 
to a recent conjecture the relative ageing of the twins is determined neither by their 
velocities (with respect to a static observer) nor their accelerations (in the case of non- 
geodesic motions). The twin travelling on a radial geodesic in Schwarzschild spacetime 
is at the reunion the oldest one among all the twins. Yet the twin moving on a geodesic 
circular orbit is younger than the static twin. 

Acknowledgements. I am deeply indebted to Zdzislaw Golda for solving a differ- 
ential equation with the aid of Mathematica and to Sebastian Szybka for making some 
analytic and numerical computations. This work was supported by a grant from the 
John Templeton Foundation. 
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Appendix 

Here we solve equations (45) and (47) and get a generic Jacobi field on the timelike cir- 
cular geodesic B what allows one for the search for all conjugate points on this geodesic. 
To solve for Z\ and Z 3 one first introduces a dimensionless independent variable x = qs 
(q is defined after eq. (43)) and then replaces the two coupled equations (45) and (47) 
by one equation for Z\ by eliminating Z 3 . As a result one gets a fourth order equation 

Z[ IV) 4 cot 2x Z"' + ( —= (3 -2] Z" + — \8 + (8- 2) cos 2ar] Z x ' 

\ sin 2x J sin 2a; 

4 
-2^- {(3 cos 2x + (3 - 2) + 3 - (3 



Zi = (A.l) 

.sin 2x 

with Z 1 = dZi/dx etc. The equation has four linearly independent special solutions^ 

Zin, N = 1, 2, 3, 4 and these read 

Z11 = sin a;, (A. 2) 

Z u = 2cosx + f3xsmx, (A. 3) 



2 sin a; sin (a/4 - (3x) + a/4 - (3 cos a; cos (a/4 - ]3x), 3 < (3 < 4, 
Z\z — { x cos x + x 2 sin x, 13 = 4:, (A. 4) 

2 sin a; sinh(A//3 - ~ix) + y/8 - 4 cos a; cos1i(a//3 -4a:), /3 > 4, 

2 sin a; cos(v / 4 — (3x) — a/4 — ]3 cos a; sin (a/4 — /3 a;), 3 < [3 < 4, 
, u -^ 4a; 3 sin a; + (3 + 6a; 2 ) cos a;, /3 = 4, (A. 5) 

2 sin a; cosh(v//3 - 4x) + a//3 - 4 cos a; sinh( A /5^ r 4x), (3 > 4, 

x = gs. One sees that (3 = 4, corresponding to r$ = 6M, the innermost stable circular 
orbit (ISCO), is distinguished. The general solution to (A.l) is then 

4 
Z\= 2_^ Cn Z\n 

N=l 

with arbitrary constants Cn- The scalar Z 3 is determined from eq. (45) as action of a 
linear operator on Z\, 

( 9 J 1 -2 \ 

cot 2x - . Zi = LZx, (A.6) 



/3 sin 2a; (ia; 2 /3sin2x 



or Z 3 = J2 C N LZ 1N = J2 c nZ 3N . (A. 7) 

N=l N=l 



5 The four solutions were found by Z. Golda by manipulating with the equation with the aid of 
Mathematica. There is a first integral generated by the timelike Killing vector field (we discuss such 
integrals in Part II) but it turns out to be of little use. 

15 



Explicitly, the special solutions Z 3N are 

Z 31 = — cos a;, (A. 8) 

Z 32 = 2 smx — f3x cos x, (A. 9) 



A33 



-2 cos a; sin (a/4 - px) + y/4 — f3 sin a; cos (a/4 - p*x), 3 < /3 < 4, 

a; sin a; — a; 2 cos a; , /3 = 4, (A. 10) 

a//9 — 4 sin a; cosh(y//3 — 4rr) — 2 cos a; sinh(\//3 — 4a;), /3 > 4, 



-2 cos a; cos(a/4 - px) - V 4 - p sin a; sin ( a/4 - px), 3 < /3 < 4, 
: :u = ^ 3(1 + 2x 2 ) sinx-4x 3 cos a;, fi = 4, (A.ll) 

\//3 — 4 sin a; sml^A/p 1 — 4x) — 2 cos a; cosh( v / /3 — 4x), (5 > A. 

Altogether the general Jacobi field on B depends on 6 arbitrary constants and takes the 
form 

4 

Z^s) = J2 C N (Z 1N e t £ + Z m e 3 l ) + (C>smp s + C"cosp s)e 2 l . (A.12) 

N=l 

We have seen in section 3.2 that the special Jacobi field proportional to e 2 generates 
an infinite sequence of points conjugate to P of which the first two are located on the 
segment PoP\, hence the field does account for the fact that the circular geodesic B is 
shorter between P an d -Pi than the radial geodesic C and the non-geodesic curve A. 
We now establish whether a Jacobi field spanned on the basis vectors el and e 3 may 
also give rise to a conjugate point. The procedure is simple and laborious. One must 
study separate cases 3<f3<A, (3 = A and j3 > 4. In the special case j3 — 4 i.e. the 
ISCO geodesic line, the initial condition -Zi(O) = = Z 3 (0), yields C 2 = 3Ci/4 and 
C4 = — Ci/2 with arbitrary C\ and C 3 . It turns out that the scalar fields Zi(qs) and 
Z 3 (qs) have no common roots for s > 0. In the generic case, /3 7^ 4, one again gets that 
the functions Z\ and Z 3 vanishing at s = depend on arbitrary C\ and C3, now with 

C 2 = ~\0-4\ 1 '*C 3 , C, = - l -C x . (A.13) 

For /3 > 4, i.e. for 3M < ro < 6M, one again finds that Zi and Z3 do not possess 
common roots for s > 0. Finally, for r > 6M one finds that the scalars 

Z x = 2Z U - Z u and Z 3 = 2Z 31 - Z u (A. 14) 

vanish at 

2 ™Mt7 1^77 . n = l,2,...,. (A.15) 



q^/A^fi \M Vr -6M 
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Accordingly, besides the sequence (50) found previously, there is a second infinite se- 
quence of points Q' n (s n ) conjugate to P on the geodesic B. To check whether the first 
of these, Q[(si), is located on the segment PqP\ we compare the value of s\ with the 
length sb of the curve from Pq to P\. From (8) one gets 



S\ ( r \ : "- 



s B \r - 6M 



> 1 (A.16) 



showing that all the points Q' n lie beyond the relevant segment. One concludes that 
these conjugate points do not account for the fact that the circular geodesic B has a 
non-maximal length. 
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